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Abstract. Conventional perturbative QCD calculations on the production of a heavy quark 
“iT” consist of two contrasting approaches: the usual QCD parton formalism uses the zero-mass 
approximation (rnn = 0) once above threshold, and treats H just like the other light partons; 
on the other hand, most recent “NLO” heavy quark calculations treat mu as a large parameter 
and always consider H as a heavy particle, never as a parton, irrespective of the energy scale of 
the physical process. By their very nature, both these approaches are limited in their regions of 
applicability. This dichotomy can be resolved in a unified general-mass variable-flavor-number 
scheme, which retains the mn dependence at all energies, and which naturally reduces to the two 
conventional approaches in their respective region of validity. Recent applications to lepto- and 
hadro-production of heavy quarks are briefly summarized. 


1 Introduction 

The production of heavy quarks in photo-, lepto-, and hadro-production processes has be¬ 
come an increasingly important subject of study both theoretically and experimentally. For 
a comprehensive review and references, see Ref. [1]. The theory of heavy quark production in 
perturbative Quantum Chromodynamics (PQCD) is considerably more subtle than that of 
light parton (jet) production because of the additional scale introduced by the quark mass. 
Let us consider the production of a generic heavy quark, denoted by H, with non-zero mass 
rriH, in high energy interactions. For definiteness and simplicity, unless otherwise stated, 
we shall use deep inelastic lepton-hadron scattering as the talking example. A reasonable 
criterion for a quark to be called “heavy” is mn ^ tigcD, so that perturbative QCD is 
applicable at the scale rriH- Thus, conventionally, {c,b,t} quarks are regarded as heavy. 

The relevant energy scales of this problem are: (i) a typical small scale such as Aqcd 
or masses of light mesons, nucleons, ...; (ii) the highest energy scale E or a/s; (hi) a typical 
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large scale in the physical process, such as pt of the heavy quark (or the associated heavy 
flavor hadron), Q of deep inelastic scattering or Drell-Yan processes, or some large mass (such 
as mw,z,Higgs, susy) - to be denoted henceforth collectively as Q] and (iv) the heavy quark 
mass niH- By dehnition, uih 3> Aqcd] and we need ^/s to be fairly large compared to uih for 
the production cross-section to be substantial. Thus the important ratio of scales remaining 
which determines the physics of the heavy quark production process is that between mn and 
Q} We shall be mainly concerned with c and b quarks for which this ratio can vary over a 
wide range in practice. 

2 Two Contrasting Conventional Approaches 

The two conventional approaches to heavy quark production in PQCD can be summarized 
by the following contrasting master equations used in the calculation^ 

ZMVFN: aiA^cx= E flCa, p) ® hia^cxC,Q, p) (1) 

a = all active partons 

FFN: aiA^HX= E flCa, p) ® ^la^nxC^Q^mn, g) (2) 

a = light partons only 

The Zero-mass Variable-flavor-number (ZM-VFN) scheme formula, Eq.l, is used 
routinely in most high energy calculations: in global analyses of parton distributions, from 
EHLQ [2] to MRS [3] and CTEQ [4], as well as in all analytic or Monte Carlo programs 
for generating SM and new physics cross-sections. In this equation, the parton label a is 
summed over all possible active parton species] g is the factorization and renormalization 
scale; and aia^cx is the perturbatively calculable hard cross section involving partons only. 
“Active” partons include all quanta which can participate effectively in the dynamics at the 
relevant energy scale p (~ Q) [2,5,6], ineluding charm and bottom quarks at current collider 
energies. Thus, the active flavor number ng depends on the energy scale (“resolving power”) 
of the problem; it is not hxed at any particular value. The hard cross-section aia^cxC, Q, g) 
is calculated in the limit of zero mass for all the partons, and it is made infra-red safe by 
dimensional regularization in the MS scheme - hence the name Zero-mass Variable-flavor- 
number (ZM-VFN) scheme. 

The advantage of the ZM-VFN scheme is that it is quite simple to implement. For the 
light partons a ={g,u,d, s}, ma —0 is a valid approximation for all hard scale Q (since, 
by dehnition, Q 3> roa)- But for a heavy quark H, it is a reasonable approximation only in 
the high energy regime g ^ Q ^ hih] and it clearly becomes unreliable in the intermediate 
region Q Oimn). 

In this talk, we shall not consider “small-x” problems associated with logarithms of the large ratio Vs/Q, 
cf [12]. 

If a final state particle C is observed, the factorization formula should also contain a fragmentation 
function dC{z,g). We leave ont d^{z,g,) here only for simplicity of discussion. All statements concerning 
the parton distributions also apply to the fragmentation functions, if present. 



In contrast to the above, the fixed-flavor-number (FFN) scheme, Eq.2 has been used 
in most recent fixed-order perturbative calculations of heavy quark production [7-9]. In this 
scheme, by definition, only light partons (e.g. m, d, s and g for charm production) are included 
in the initial state: the number of parton flavors Uf is kept at a fixed value regardless of the en¬ 
ergy scales involved (n/ = 3,4 for c, b production respectively). The main feature here is: H 
is pictured as a heavy particle - much in the same way as W, Z, and other new heavy particles, 
and very different from the zero-mass light partons - hence the mass rriH is kept exactly in 
the hard cross-section Q) ^h, h)- Typically, the perturbative Qi ^h, h) 

will contain logarithm factors of the form a^{p) /niH) If Q ~ mu (and 

X ~ Q'^/s is not too small), these factors are under control; and we have effectively a one 
large scale hard process. Hence, the FFN scheme is the natural scheme to use in the energy 
region Q ~ niH - this is precisely where the ZM VFN scheme is expected to be inappropriate. 

From the heavy particle perspective, this approach also has the advantage of being concep¬ 
tually simple, even if the NLO calculation requires considerable amount of work. However, 
the sharp distinction drawn between the H quark and the other light quarks, say between 
c and s, in this formalism appears quite unnatural as the hadron system is probed at the 
scale /r ~ Q available in current high energy processes. And it has been known since the 
next-to-leading order (NLO) calculations in the FFN scheme were completed [8,9] that, for 
both charm and bottom production, there are two disconcerting features about the results: 
(i) the NLO corrections turn out to be of the same numerical magnitude as (in fact, gen¬ 
erally larger than) the leading order (LO) result; and (ii) the uncertainty of the theoretical 
calculation, as measured by the dependence of the calculated cross section on the unphysical 
scale parameter p, is as large in NLO as in LO - contrary to what is expected from a good 
perturbation expansion [10]. These features mean that the truncated perturbative series in 
this scheme has left out important physics effects. Experimentally, it is also known that the 
measured charm and bottom production cross sections do not agree with the NLO theoret¬ 
ical predictions, at least in the overall normalization, even when the scale fi is allowed to 
vary within a reasonable range [11]. 

This situation may not be all that surprising: for c and b quarks, the condition Q ~ mn 
is not well satisfied in most practical cases. In fact, current experimental ranges for lepto- 
and hadro-production of these heavy flavors mostly lie in a region between those appropriate 
for the ZM VFN ( Q 3> mu) and FFN {Q ~ mH) schemes. We need a well-defined theory 
which applies over the full Q range! Other possible sources for these problems are: (i) large 
corrections due to large logarithms of (s/Q^)—the small-x problem [12]; (ii) inadequate 
understanding of the hadronization of heavy quarks in comparing PQCD calculations with 
experiment; and (iii) existence of non-perturbative components of H inside the nucleon which 
are, by definition, excluded by the FFN scheme. In this talk, we shall concentrate on physics 
issues pertaining to the changing role of the heavy quark H over the full Q range. It is 
particularly interesting because the interplay between the two independent scales mn and 
Q embodies much interesting QCD physics which is amenable to precise treatment. 



3 A Unified, General-mass, Parton Approaeh 


When the energy scale becomes large, Q/mn 3> 1, the FFN scheme becomes suspect 
because large logarithm factors of \n{Q/mu) in the hard cross section becomes in¬ 

creasingly singular, and higher-order terms containing higher powers of the same can no 
longer be omitted. In other words, the truncated perturbation series in this scheme can be¬ 
come rather unreliable as Q becomes large. The clue for addressing this problem is already 
contained in Eq.l: large logarithms of the form a” \n"'~^{Q/mH) in dfiGfix can be resummed 
to all orders in to become the parton distribution (evolved to {k + l)-loops). 

The H parton should be included in the sum over parton flavors; it participates in the hard 
scattering on the same footing as the other partons. After removing these potentially dan¬ 
gerous logarithm terms, the remaining hard cross section dab^HX becomes infra-red safe as 
Q/niH —^ oo. It is important to note, however, the resummation of large \n{Q/mH) loga¬ 
rithms does not require taking the mn —^ 0 limit for the remaining (“mass-subtracted”) hard 
cross-section as is done in the conventional ZM VFN formulation, Eq.l. In fact, by retaining 
the rriH dependence in the mass-subtracted (hence infra-red safe) dab^uxU^Q^mH, p), one 
arrives at a consistent theory for heavy quark production which is valid over the entire energy 
range from Q < mn to Q 3> mH3 

aiA^cx{s,Q,mH) = fAGa,p)®dia^cx{s,Q,mH,p) (3) 

a = all active partons 

A program to systematically implement this intuitive physical picture has been developed in 
a series of papers in [6,13-15]. The resulting formalism constitutes a natural generalization 
of the conventional zero-mass QCD parton formalism to correctly include general quark mass 
effects, hence will be called the general-mass variable-flavor-number (GM-VFN) scheme. (In 
some recent literature it has also been called the ACOT scheme. Ref. [15].) 

More precisely, this formalism is based on a well defined renormalization scheme [16] which 
provides a natural transition from the threshold region Q ~ 0{mH) to the high energy region 
Q ^ rriH] and the validity of the generalized factorization theorem can be established order- 
by-order in perturbation theory [17]. The key points are: 

• the renormalization scheme is a composite of two simple schemes, natural for Q "A oih 
and Q > mu respectively, with matching conditions that make the schemes equivalent 
in the domain of overlap Q ~ mu where they are equally valid for practical low order 
calculations [16]; 

• one scheme utilizes a subtraction procedure (BPHZ) which leads to manifest decoupling 
of the heavy particle in the region Q ttih, thereby gives precise meaning to the FFN 
scheme (with no heavy quark partons); 

• the other scheme is ordinary MS as regards the definition of the coupling as{p) and 
the parton densities f'^{x,p), hence retains the normal {mu = 0) evolution equations 

Factorization of any applicable fragmentation functions is implicitly assumed. 



for the latter in the region Q > tuh [6] - this conies about because the evolution 
kernels are anomalous dimensions which are derivatives of renormalization constants, 
and renormalization constants are mass-independent in the MS scheme; 

• the factorization scheme is dehned such that all infra-red safe m/f-dependent effects 
are preserved in the hard cross-sections, so that there is no loss of accuracy when 
Q ~ mn [15]. This is accomplished by dehning d'ia^cx{s,Q,mH, as the full 
<Jia^cx{s,Q,mH, with mass {mn) singularities subtracted. 

These features guarantee that predictions of this formalism: (a) coincide with those of the 
FFN scheme in its region of applicability, Q ^ uih', (b) reduce to those of the conventional 
zero-mass parton model in the asymptotic energy regime Q 3> uih] and (c) provide a good 
approximation to the physical cross-section over the entire energy range in between, since the 
remainder of the perturbation series contains no large logarithms of the kind log{Q/mH)-^ 


4 How, and Why, does the ACOT scheme work? 


How and why does this scheme work were described in some detail for lepto-production 
of heavy quarks in Ref. [15]. We recall the essential points here. Writing out the hrst two 
terms in the perturbative expansion of Eq.3, we have 

aiA^Hxis,Q,mH) = fA ® + Ia® ( 4 ) 


where the superscript n in denotes the order in Og of the hard cross-section a. The lowest 
order hard cross-section is identical to the Born expression since the tree 

diagram does not need any subtraction. The order Og hard cross-section is given by 
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where ^ is the perturbative parton distribution function of Ending H in g] PgH is the 
usual g ^ H splitting function; and ^crig^i^ff is the unsubtracted cross-section^ for the 
gluon fusion process Ig IHH. The subtraction term above can be formally derived by 
applying Eq.4 to the partonic cross-section ^(Jig^i^jH (with light-parton colinear singularities 
subtracted), and then solving for [15]. 

Substituting Eq.5 in Eq.4, the physical cross-section becomes: 


aiA^HX — f A ® + f A ® ^ f 4^ ® 

( 6 ) 

= Ua - Ja ®^ /f ) ® + Ja® 

y Needless to say, colinear singularities associated with light partons are MS subtracted. 

As mentioned earlier, we do not consider “small-a;” corrections in this paper. 

Here “unsubtracted” refers to heavy quark mass effects. As mentioned before, colinear singularities due 
to light partons are always subtracted as in the MS scheme. 



The three terms on the right-hand-side of the hrst line are: heavy-flavor excitation (HE), 
heavy-flavor creation (HC), and the subtraction term; the last physically represents the over¬ 
lap between the two production mechanisms. Both and ^(Tig^iuR (unsubtracted) 

contain the full mn dependence. For Q 3> uih, it is useful to view the right-hand side 
of Eq.6 as on the hrst line. The quantity in the parenthesis (really just is free 

of large In logarithms because of the subtraction; it is infra-red safe. The whole term 
remains at numerical order Ug x 0{1) (in contrast to the unsubtracted ^(Xig^iuH which is the 
LO FEN scheme result with a large log factor, ag x (9(ln^-)) in the large Q limit. Thus, 
one recovers the LO ZM VFS formula with the HE contribution as the dominant term for 
the cross-section. In the threshold region, Q ~ uih, it is more useful to focus attention 
on the second line in Eq.6 : the two terms inside the parenthesis are individually small in 
this region and, in addition, they cancel each other up to order since both satisfy the 
evolution equation (to this order) with the same boundary conditions (assuming there is 
no non-perturbative charm). As a consequence, the full cross-section is dominated by the 
second (HC) term - which is the FEN scheme result to order ag. 

The explicit expression for the subtraction term clearly shows how it overlaps the HE and 
HC mechanisms, hence leads to the appropriate cancellations in the respective kinematic 
regions. The same principle applies in other lepton-hadron and hadron-hadron processes to 
all orders of perturbation theory [17]. Numerical calculations based on Eq.6 conhrm the 
features described above for relevant physical cross-sections. See Ref. [15] and the talk by 
Schmidt [18] for plots of {x, Q) with individual contributions from the three terms on the 
right-hand side of Eq.6 which explicitly illustrate these features. 

5 Complementarity between HE and High-order HC Contributions 

The GM-VFN (ACOT) scheme highlights some overlapping features of HE and higher or¬ 
der HC mechanisms for heavy quark production (hence the need for the subtraction to avoid 
double counting): the two are not mutually exclusive, as sometimes perceived; rather, they 
are complementary. For the total inclusive cross-section (e.g. structure functions {x, Q) 
in DIS), in particular, the HE contribution represents the result of resumming the colinear 
parts of HC diagrams to all orders in the running coupling. This provides an efficient method 
of obtaining important quantitative results without having to calculate many complicated 
higher order diagrams in the FEN scheme. 

The trade-off is that some information on the differential distributions (e.g. transverse 
momentum spectrum of the heavy particle) is integrated over in the resummation, hence 
becomes less accessible. For example, if one is interested in the pt distribution of the charm 
quark in DIS lepto-production with respect to the virtual photon-target axis, the HE con¬ 
tribution, as well as the subtraction term, in Eq.6 will be formally proportional to the delta 
function 5{pt). This is, mathematically, a distribution (rather than an ordinary function) 
which needs to be folded with a suitable smearing function - some combination of experi¬ 
mental resolution function and theoretical pt distribution (see below) - in order to produce 
meaningful physical cross-sections. Away from the small pt region, the low-order FEN scheme 



diagrams will be the logical place to start in obtaining the leading contributions to the phys¬ 
ical pt spectrum. Cf. the talk by Schmidt [18]. To obtain a precise theory of pt distributions 
over the entire range, a different kind of resummation will be required. Thus, for heavy 
quark production in PQCD, as in many quantum mechanical problems in general (say, the 
double-slit problem), the appropriate way to formulate the theory depends intimately on the 
physics question asked. 


6 Recent Developments 

The ACOT scheme has been applied to the analysis of recent charm production data 
in neutrino scattering [19], yielding the most up-to-date information on the strange quark 
content of the nucleon. 

Recently, in the wake of new precision measurements of the total F 2 [x, Q) at small x [20] 
(where charm hnal states consist of up to 25% of the cross-section) as well as first mea¬ 
surements of F 2 {x, Q) [21], it has become obvious that a more precise theoretical treatment 
of charm production in NC DIS than those used so far (the ZM-VFN and FFN schemes 
described earlier) is now needed. A hrst step is taken by the CTEQ group, performing a new 
global QCD analysis of parton distributions [22] based on the GM-VFN scheme of ACOT 
(to order only). A similar study has been done by MRRS [23], using a related procedure 
(see also [24]). An order calculation in the general mass scheme is within reach [25] since 
much of the known FFN scheme results to this order can be used as the starting point. 

Of equal importance is an improved theoretical treatment of various heavy quark cross- 
sections in hadron-hadron scattering beyond the conventional calculational schemes, partic¬ 
ularly the NLO FFN scheme results which suffer from the problems mentioned in Sec. 2. 
Some preliminary results have been obtained in the ACOT scheme [26,27]. The qualitative 
features of these results, compared to existing calculations, are similar to those on lepto- 
production described briefly in Sec. 4. As expected, the inclusion of heavy-flavor-excitation 
contributions leads to: (i) reduced dependence on the (unphysical) scale parameter /r; and (ii) 
an increase in the predicted cross-section over most of the range of pt of the produced heavy 
quark, as preferred by data [11]. Detailed phenomenology of the inclusive cross-section, as 
well as correlations of hnal-state heavy particles still need to be pursued [27]. 

7 Concluding Remarks 

The general-mass variable-flavor-number scheme of Ref. [6,14,15,17] generalizes the famil¬ 
iar zero-mass variable-flavor-number QCD parton framework (valid only at large Q scales) 
to include quark mass effects, and it reproduces the results of the widely used FFN scheme 
heavy quark calculations (valid in the one-large-scale region Q ~ m/f). It naturally unihes 
the two contrasting conventional approaches in a well-dehned renormalization and factoriza¬ 
tion scheme. First results on lepto- and hadro-production demonstrate improvements over 
existing calculations both in smaller /^-dependence and in increased cross-sections. 



However, the more complete theory is not a cure-for-all. Although detailed phenomenology 
has yet to be done, there is no doubt that further developments in several directions are 
needed to reach a full understanding of heavy quark production in QCD. For instance, (i) 
the general mass scheme should be implemented to higher orders for both lepto- and hadro- 
production; (ii) the signihcance of possible large logarithms of the type ln((5^/s), ln(m|^/s) 
~ Inx - the small-x problem - needs to be better understood [12]; (iii) both perturbative and 
non-perturbative aspects of the hadronization of the heavy quarks deserve further study; and 
(iv) the question “Are there non-perturbative charm/bottom components inside hadrons?” 
needs to be answered. The last question, having been discussed in the literature for many 
years, has to be carefully investigated phenomenologically; and this can be done only in the 
framework of the general-mass scheme, since the existence of non-perturbative non-zero-mass 
parton is excluded by assumption in the FFN scheme. 
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